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Analysis of state reconstruction both classical and quantum mechanical on equal footing is out- 
lined. The meaning of "mutual unbiased bases" (MUB) of Hilbert spaces is explained in detail. An 
alternative quantum state reconstruction, that utilizes mutual unbiased bases (MUB), is given. The 
MUB approach is then used for state reconstruction in a finite, d, dimensional Hilbert spaces. 
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\ The position operator, x and the momentum operator p form a complete operator basis [l|, i.e. an operator that 
commutes with both must be a scalar. In view of this one might expect that the reply to W. Pauli query 0, viz: can 
\ one construct a wave function, ip(x), amplitude and phase, from the probabilities for (all) x and p, would be in the 
affirmative, i.e. that ^(x)! 2 and |V>(p)| 2 have sufficient information to deduce ^{x). However this is not sufficient. 
' Indeed what is required to reconstruct a state is the issue dealt with in this pedagogical article. In other words we 
will see that measuring (assuming unlimited ensemble of |^)) the position (x\tj)) (ij)\x) and the momentum (p\?p}(ip\p) 
would not allow the reconstruction of the state. To reconstruct the state one requires the probability distributions of 
all the so called mutually unbiased operators (these will be defined below) - this set, of which the operators x and p, 
are members will now be considered. 

The conceptual idea of mutual unbiased operators was introduced by Schwinger [l[ when he considered vectorial 
bases for Hilbert spaces which exhibit "maximal degree of incompatibility" - the eigenvectors of x - \x) and of p, \p) 
are example of such bases. The information theoretical oriented appellation "mutual unbiased bases" (MUB) were 
O" 1 introduced by [H, IE 0] with the following formal definition: Two orthonormal vector bases, Si, 2?2, are said to be 
mutually unbiased (MUB) if and only if (B\ ^ B2) 

j> ■ V |m2) e 81, Bi resp. \(ui\iL2)\ = constant, (1-1) 

£f) i.e. the absolute value of the scalar product of vectors from different bases is independent of the vectorial labels within 
Tj" ■ either basis. This implies that if the state vector is measured to be in one of the states e.g. \ui) of the base Bi, it is 
I equally likely to be in any of the states \u2) of the base £> 2 . (The value of |(iti|«2)| ma Y depend on the bases, B\, B2 
. as indeed is the case for the continuous dimensionality - d — ► 00, - limit.) MUB were found of interest in several 
7— I ' fields. Thus the ideas are useful in variety of cryptographic protocols (e.g. pj), signal analysis 0. They are of partic- 
CTN , ular interest in quantum state tomography, Q where quadrature fl"ot | observations allow the construction of Wigner 
■ function of the state under study. The relation of MUB to Wigner function is extensively studied by Wootters and 
^! \ co-workers [E [HQ. Much of these works were devoted to finite dimensional Hilbert spaces. Here it was shown that a 
V~j ■ d dimensional Hilbert space can accommodate at most d+1 MUB and the value of |(ui|u2)| is -j= [1, IbH [TR . 
The finite dimensional theory is intriguingly connected to algebraic field theory [E [E EE EE EMEE HI ; connections 
to other sophisticated mathematical notions are given in [2l|, [22|, HE HE HH and summarized in [l7|. Extensive infinite 
dimensional studies are given in 0, [E [E Hfij . We will argue that the information required to reconstruct a quantum 
state may be gained by measuring all the mutually unbiased operators. 

The problem of state reconstruction has, of course, its classical analogue. Within classical physics we may refer to the 
density in phase space, p(x,p), of a physical system as the state (of the system). In the following we will study how, 
i.e. what kind of measurements, will suffice to reconstruct such a state. The methodology will be seen to based on 
the inverse Radon transform Q. Then, utilizing the quasi probability nature of the Wigner representative function 
W p (x,p) [E Q of the quantum state operator p, [ll[, we show that reconstructing the Wigner function is closely 
analogous to the reconstruction of classical phase space distribution function. This is the content of Section II. 
Section III contains an alternative derivation of the quantum state reconstruction which is based on MUB. This 
approach does not require an explicit use of Radon transform. In Section IV we give an MUB based state reconstruc- 
tion for finite, d, dimensional cases with d a prime number. This approach does not involve reconstruction of the 
Wigner representative function whose finite dimensional version is rather challenging [J, EH ■ The last section contains 
summary and discussion. 
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II. CLASSICAL-QUANTUM STATE RECONSTRUCTIONS ANALOGY 

The classical state, p(x,p), is the phase space distribution function. In terms of this state the classical probability 
(density - to be understood henceforth), p{x' , p'), for the particle position to be at the local x = x' and the momentum 
of p = p' is 



p(x',p')= J dxdpp{x,p)5(x — x')S(p — p'), (2-1) 
The Fourier transform of the distribution, F(a, b; p), is 

F(a,b;p)= J ^ e ™*e ib *>p(x,p). (2.2) 



Thus, 



2k 

The marginal probability for the particle being at x = x' while being completely ignorant of its momenta, is 



p(x',p') = I < ^e- iax 'e- ibp 'F(a,b;p)= I dxdpp(x,p)S(x - x')6(p - p 1 ). (2.3) 



p(x',0)= dxdpp{x,p)5{x~ x'). (2.4) 



(The entry label will become meaningful below.) Correspondingly the marginal distribution that the particle's phase 
space coordinate be x' = Xg = xC + pS; (C = cosQ, S = sinO) is given by the following transform, 

p(x',6)= [ dxdpp(x,p)S(x' -xC-pS). (2.5) 



This transformation is known as Radon transformation, Q. p(x' ,6) is a directly measurable quantity 0]. Before 
considering the inversion of this, i.e. getting p(x,p) in terms of p(x',6), we note that such, (Eq.f |2~5"|) ). marginal 
distribution implies complete lack of knowledge of the odds for the realization of phase space coordinate p' = Pg = 
—xS + pC . The proof is as follows: define the 9 depending variables 

Xg = Cx + Sp, pg = —Sx + Cp. (2-6) 

Now one readily verifies that 

p(x',6)= / dxdpp(x,p)S(x' - xC - pS) — dxgdpgp(x(xg,pg),p(x e ,p e ))8(x' - xg). QED (2.7) 



The inversion of the formula proceeds as follows. Expressing the delta function in Eq.f 12.5) ) in terms of its Fourier 
transform 



p(x',6) = J dre-* xr F(r;6,p), 
F(r;6,p) = [ dxdpp(x,p)e tCrx e lSrp , i.e. 



F(r;0,p) = F(a = rC,b = rS; p), and, 

Hr;6,p) = J ^- e ix ' r p(x' ,9). (2.8) 
This relates the measurable p(x' , 9) to F(r; 9, p). The inverse Radon transform is obtained as follows 

p{XjP) = [ ^ e -^ e -^ F{aMp) = r rdr r ^-e- irC *e- irS *>F(r;9,p). (2.9) 
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Performing the r integration and noting that 

p(-x',e)=p(x',e + ir) (2.10) 

we get 

p{x > p) = -^l L (xc +P s-*r (2 ' n) 



This formula reconstructs, via the inverse Radon transform, the full (classical) distribution, p(x,p) from marginals, 
"sliced" distributions, p(x',8), which may be obtained by direct measurements. As noted above, Eq. (|2.5p is an 
example of (the direct) Radon transform Q. 

The quantum reconstruction utilizes the close analo gy between the classical phase space distribution function and 
the Wigner representative function W p (x,p), (3.l9Ll27Ll29ll30j| of the density operator p, 

W p (x,p) = J ^y{x- V -\p\x+ V -). (2.12) 

Consider the operator Xg = xC + pS and its eigenstates, 

X e \x',6) = x'\x',6), (2.13) 

these form a complete orthonormal basis [J]. Now the probability of finding the particle in the state \x',6) given that 
its (the particle's) state is p is, 

p Q (x',6) = Tr\x',6)(x',6\p 

= {x',6\p\x',e). (2.14) 

The quantity (x' ,9\p\x' ,9) is directly measurable Q. Transcribing this to the Wigner representative function, 

P Q (x',9) = J ^W w , e) {x,p)W p {x,p), (2.15) 
Noting that (the proof is given in the next section) 

, „ „ 1 



We evaluate, 



-.e 



^([x 2 +x' 2 ]C-2xx') ^ (2.16) 



Eq. (|2.14|) now reads, 



I ± ye wy {x ^ V_\\ x ^ e){x ^ e \ x+ y_ ) = 5 {x' -xC-pS). (2.17) 



/ ^k S (x' -xC- P S)W p (x,p). (2.18) 



Which, using the analysis given above for the classical case, Eq. (|2.8l2.9[) . gives, @ 



wl , v r f +oc p Q W,e)dx'd9 

Wp[x > p) = -^Jo L (xc+ps-x'r (2 - 19) 



Thus measuring (x' ,9\p\x' ,8) i.e. pq(x,9) (cf. Eq. (|2.14p . yields the Wigner representative function for the state p 
thereby reconstructing the state since W p (x,p) is a faithful representative of the quantum state @, H, 
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III. STATE RECONSTRUCTION VIA MUTUAL UNBIASED BASES (MUB) 

We now outline an alternative state reconstruction strategy which is based on MUB. We begin with a brief review of 
a recently [1,[28[ MUB construction: The complete orthonormal basis {\x', 9)} as defined above, Eq. (|2.13[) . constitute 
a set of bases labelled by 6, each complete and orthonormal.. We shall now argue that these arc MUB. Wc have, as 
can be checked directly [2, [28| , 

X e = U*(0)xU(0), (3.1) 

with, 

a = ±(x + ip); ^ = -L(x - ip); U{9) = eT^\ (3.2) 

Returning to Eq. (|2.13p . wc have that solutions of the equation may be written in terms of {|a;)}, the eigenvectors 
of the position operator, as 

\x; 9) = U\6)\x). (3.3) 
This defines our phase choice. Thus we may read off the x-representative solutions [28l. [3lj. 

,([x 2 +x' 2 ]C-2xx') ^ /g,4) 



( x '\ x -8) - 1 e -A( u2 ^' 21 
( 1 ' } ^\S\ 



This can be verified directly as the solution of 

(Cx-iS^-){x'\x,e) =x(x'\x,6). (3.5) 
v ox 

With the phase choice above the state, x «-> x', is symmetric for property that facilitate several 

calculations. We can now verify that the bases {\x; 9)} and {\x / ; 9')} with 9^9' are MUB: thus [3, HI] 

(x'-9'\x-9) = (x'\U\9-9')x)\ = (x\x;6-6') = , 1 e -^hPi^ 2 +^ c ^fi')-^') , ( 3 . 6 ) 

^2k\S(9,9')\ 

Here S(6, 6') = sm(9 - 6') and C(6, 9') = cos{9 - 9'). Indeed, 

\{x'-9'\x-9) \ = - 1 (3.7) 

the number | (x'\ 9'\x; 9) \ is independent of the vectorial labels x and x' , i.e. the bases {|a;;0)} with distinct 9 are 
MUB. We further note that 

\im(x\x';9) -> S(x - x') 

}im(x\x';e) - ^L. (3.8) 

This concurs with the observation that, (c.f., Eq. (2.4)), in the first limit |x;0) is the eigenfunction of x whose x 
representative is 8(x — x'). In the second limit, \x'; is the eigenfunction of the momentum operator (p) thus the 
limit is simply a Fourier representative expressing the momentum in the x representation. 

We now derive an explicit tomographic formula for the density operator i.e. the state is reconstructed via expressing 
the state, (x\p\x') in terms of measurable (x, 9\p\x, 9). To this end we note that the set 

e M V bj5 , -oo<a,6<+oo, (3.9) 

is a complete and orthogonal operator basis fl3j 

dadb 



2- 



■tr(e < °*e a *[e <a '*e tt '*] t ) = 5{a - a')S(b - b'). (3.10) 
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Thus wc may expand, 



— / dadb{trpe lai e tb P)(e lai: e lbp y. (3.11) 

27T J 



Defining a = rcos9, b = rsinO and using the equality. 



giax^ibp g(i r s ^ 2 ) ^i(cos6x+sin&p)r 



we may rewrite Eq. p. lip as, 

p - J°° J* rdrd0 ( trpc (i£M^) c i(cosei&+8in6ii)rj ^{i^^) e i{cosex+sinep)ry _ 



10 J — 7T 

Utilizing the spectral representation, 



i(cosex+sinep)r = dc\c, 9}e lrc (c, 9\, (3.13) 



and evaluating the trace via the orthonormal bases {|c; 6)}, this expression becomes 



00 rdrd9 

dede' I / — {c;e\ P \c;9)e %rc \c';e)(c';6\e- tcr . (3.14) 

J-7T 27r 



Evaluating the integral 

/>oo J poo -t 



(3.15) 



and inserting it in the above expression we get our tomographic formula (the bounds on the 9 integration is justified 
via Eq.flHTO])): 



P = — dede' / -^—( C ;e\p\c;e)\c';6)(c';6\. (3.16) 
n J -00 Jo ( c ~ c ) 



Thus the density matrix, (x\p\x'), is given in terms of the observables (c;9\p\c 

-p d6 e 75{( x2 - x ' 2 ) c - 2 ( x - x ') c ') 

7o T^J? {cMc ' e) 



oo 



x\ P \x') = --j dede' / 7 — J ^(c;6\p\c;9)- — . (3.17) 



This completes our state reconstruction via MUB: one may construct (x\p\x') in terms of the measurables (c; 0\p\c; 6). 
The expression is of course equivalent to the one given in terms of Wigner function Q : Evaluating the Wigner function 
in Eq. 16[) gives, 

Where the Wigner function is given in terms of measurements that give the probabilities (c; 0\p\c\ 9), conforming with 
1- 

IV. STATE RECONSTRUCTION IN FINITE DIMENSIONS 

The difference in the two strategies discussed above for the continuous phase space (Hilbert space dimensionality 
d — ► 00) involve, in essence, the benefit of two ways of handling the complete orthonormal operators e lax e lb P_ Thus 
when viewing them as Fourier transformation 0] we are led to the Wigner function via the Radon transformation. If, 
on the other hand, we consider the spectral representation (note a = rC, b = rS, c = cos9, S = sin9), 

el ax +l bp = e tr(£C+pS) = j d X '\x' ,6)?™' {x' ,6\, (4.1) 
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we are headed for state reconstruction via MUB. This approach will now be broadened to accommodate the finite ,d, 
dimensional Hilbert space state reconstruction 

Schwinger [l[ noted that the physics of finite dimensional, d, Hilbert space is accountable via two unitary operators, 
X, and Z. Thus if we label the d distinct states, termed the computational basis, by \n), n = 0, 1, ..d—1; \n + d) = \n), 
these operators are defined by: 

Z\n) =cj"|n); X\n) = \n+l), (4.2) 

with uj = e 27Tt / d . They form a complete set, i.e. only a multiple of the identity commutes with both X, Z. We 
shall briefly outline a method utilizing these operators (due in the main to [ill]) to construct the d + 1 MUB for a d 
dimensional Hilbert space with d being an odd prime. This review will be of help in building our sets of entangled 
states that we shall associate with these MUB. The computational basis vectors spans the Hilbert space. All operators 
in this space are expressible in terms of the d 2 Schwinger-operators [l[: 

X m Z l ; m,l = 0,1,. ..d- 1. (4.3) 

The operators X m Z l whose number is obviously d 2 form an orthogonal basis for all operators in the d dimensional 
Hilbert space, 

= 6 mjm 'Si,i'- (4-4) 

This follows from Eq. (|4.2p which implies the commutation formula 

XZ = ujZX . (4.5) 

Now, let us confine ourselves to cases wherein m, I E where is a d dimensional Galois field. In this case, we 
can relate Schwinger operators to MUB. For this aim we group Schwinger operators (|4.3|) into d + 1 sets of d — 1 
orthogonal commuting operators (which together with the identity operator form a complete operator basis to Hilbert 
space). Each set of (commuting, orthogonal) operators defines a unique vector basis in Hilbert space. All the d + 1 
sets of bases form an MUB set. 

Let us first consider the case m = in Eq. (|4.3p . Readily, the operators Z l with I = 0, ..d— 1 (1 = is, trivially, the 
identity operator) form one set of commuting and orthogonal operators. This set is diagonalized in the computational 
basis (c.f. Eq. (|4.2j) ). Next consider the case m ^ 0. In this case a unique inverse m _1 is define on Fd and thus we 
can rewrite the operators l|4.3[) as 

X m Z l = u u (XZ b ) m , (4.6) 



Tr 



X m Z l 



(x m 'z i y 



where b = l/m, v = — ^m(m — 1) and 6 = 0, 1.., d—1: m = 1, d — 1. Thus, we have associated the d(d — 1) of the 
above operators in the following manner 

x m z bm _ (XZ h ) m . (4.7) 

That is, these operators differ at most by a unimodular c number. Now, for a fixed b, we have (d— 1) orthogonal and 
commuting operators. There are d distinct such sets, each labelled by b b = 0, 1, ...d — 1 which are orthogonal 

= d 5b,b'S m ,m', m,m' ^ 0. (4.8) 

When each b-labellcd set of the d-1 orthogonal, commuting operators is supplemented with the identity operator it 
constitutes a set of d unitary, orthogonal and commuting operators. And as such, defines a vector basis for Hilbert 
space: The d (orthonormal) vectors that diagonalize these operators in the set. (We remark, in passing, that the basis 
is defined up to a choice of a phase factor which does not affect the following results.) For each (b-labcllcd) set there 
exist a a unique vectors basis. Here, we are designating the vectors that form the basis by \b;c), where the index b 
labels the basis and c that particular vector in the basis b (b, c = 0, l...d— 1). The expressions for these states in 
terms of the computational basis is , 



Tr 



(XZ b ) m ( (XZ 



\b;c) 



(4.9) 
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The eigenvalues are uj c . 

The importance of the classification of the operators in this manner arises once considering the relation between 
the different sets of the vector basis. One can readily check that these bases form MUB, i.e., 

(b;c\b;c) = 8 C ^, 

\(b',c'\b,c)\ = i b^b', (4.10) 

Thus, these d distinct sets of bases plus the computational basis (which mutually unbiased to all of these sets) form 
the maximal number of MUB, that is d+1. The proof of the last formula involves the well known Q Gaussian sums. 

Since Schwinger operators, Eq. (|4.3[) form an operator basis for d-dimensional Hilbert space, so do the set of 
operators (XZ b ) m together with Z b (b = 0, .., d— 1; m = 1, ..d— 1). Hence we may write an arbitrary density operator 
as 



p = J(EETr[ P (izr] ((xz b y»y 

, ro=l fc=0 



d 

d-l 



j2^[ P z i ](z i y . (4.ii) 



i=o / 

It is convenient to rewrite this equation upon adding and subtracting m = terms: 

-1 d-l d-l 

1 / \ 

P 



1 / d-l d-l d-l \ 

= ^ E E Tr [p( xz T] ((xz b ry i + e Tr [ P z i ] (z i f . (4.12) 

\ m=0 6=0 1=0 ) 

We now utilize the spectral representation for each of the operators, 

XZ b = EK^ m (fr>™l ->Trp(XZ b ) m = E<Mp|M}w nm . (4.13) 

m n 

Hence, 

b.n n 

Or, 

(n'\p\n") =E( n 'l fe - m )^ n l'°l 6 ' n )( 5 '"l n ") + ( n 'l' 9 l 71 ") ,5 «'^" -<5n',n"- ( 4 - 15 ) 

This is an expression for p in terms of probabilities. The numbers (b; n\p\b; n) correspond to probabilities to find 
the state \b;n) in p i.e. observable quantities. We see that, as shown in [T3, the state of the system, i.e. p, is 
reconstructed via the d+1 measurements. Each of these measurements yields the d — 1 independent probabilities 
outcomes (since the sum of the probabilities add to one). This gives (d + l)(d — 1) = d 2 — 1 numbers that determine 
the density matrix. It should be noted that this holds because the operators are non degenerate. 



V. CONCLUSIONS AND REMARKS 



An approach to classical tomography is given in close analogy with the quantum state reconstruction scheme based 
on the Wigncr representative function obtained via the Radon transform. We then reviewed an alternative route for 
state reconstruction which is based on the notion of mutual unbiased bases. The latter was shown to be applicable 
to the finite dimensional Hilbert spaces eschewing thereby the somewhat sticky issue of finite dimensional Wigner 
function. An overview of quantum ideas: measurable quantities, Wigner function and mutual unbiased bases (MUB) 
is outline. 

The analysis underscores the intriguing fact that to reconstruct a quantum state we require the probabilities of all 
the phase space plane - not merely the probabilities along the position and momentum axes as might be implied by 
(the wrong) positive reply to Pauli's query posed at the introduction. 
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